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ABSTRACT 

We revisit the polarisation induced by Faraday rotation when Cosmic Microwave 
Background photons traverse magnetised plasma. We compute the secondary B-mode 
angular power spectrum from Faraday rotation due to magnetic fields in galaxies 
and galaxy clusters with masses ranging from 10"'^^ to lO^^'^M©. We investigate its 
dependence on the electron and the magnetic field profiles. Namely, we consider both 
the /9-profile of electron density as well as an electron density distribution based on 
the Navarro-Frenk- White dark matter profile. We model the magnetic field structure 
in galaxies and clusters motivated by recent observations. We further account for its 
redshift evolution and we examine the importance of its coherence length. We find 
that the -B-mode polarisation from Faraday rotation depends on the normalisation 
parameter C/ cx Ug"^. At 30 GHz for erg — 0.8, the _B-modes from Faraday rotation 
range between 0.01 fj,K^ and 4 x 10~^ fiK^ at I ~ 10^ in the case of a maximally 
coherent fields. For smaller coherence lengths, those amplitudes are smaller and they 
peak at higher multipoles. 
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1 INTRODUCTION 



Recent cosmic microwave background (CMB) observations allowed significant progress in cosmology. The CMB temperature 
anisotropies have been measured in detail and, as a result, we know that the universe is spatially flat and that the power 
spectrum of the density fluctuations is consistent with scale- invariance i Spergel et al. 20061 ) ■ strongly supporting the inflation 



scenario. 



One of the main targets for the next generation of CMB observations is the detailed meas urement of the CMB polari- 
sation. The CMB polarisat ion is traditionally split into 2 components: i?- modes and B-modes (jKamionkowski et al.lll997b 



Zaldarriaea fc SeliakI 119971 ). Particularly, the detection of the B-mode polarisation is important for further testing the in- 



flation scenario. Gravitational waves, predicted by i nflation scenario but not detected yet, generate the pri mordial B-mode 
polarisation, while the density fluctuations do not IjKamionkowski et al.lll997al : ISeliak fc Zaldarriagal 119971 ). Therefore, the 
measurement of B-mode polarisation is an essential step toward constraining the inflationary scenario. 

H owever there are other sources of B-mode polarisation. One of them is gravitational lensing (jZaldarriaga fc SeliakI 



19971 ) . Gravitational lensing distorts the CMB polarisation fields so that a secondary B-mode polarisation is induced from the 



primary B-mode polarisation. These secondary B-modes contaminate the primary B-mode polarisation. However, signiflcant 
effort has been made to find ways how to remove the secondary B-mode polar isation due to gravitational lensing from the 
polarisation maps (see for instance Hu fc Okamoto 2002; Hirata fc Seliak 20031 ) . in order to recover the primordial B-mode 
angular power spectrum with high accuracy. 

The presence of cosmic magnetic fields, observed on various scales, is also an important source for both primary and 
secondary B-mode polarisation. The main ways of detection and measurement of magnetic fields are Zeeman splitting effect, 
synchrotron emission and Faraday rotation. Zeeman splitting provides direct observation of magnetic fields. However, the 
signal is so small with respect to Doppler broadening that it is essentially useless in the cosmological context where magnetic 
fields are week and velocities are high. Extra-galactic magnetic fields are mainly detected by the other two methods. The 
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amplitude of the fields measured in galaxies and galaxy clusters are typically of the order of 1-10 /iGauss. These methods only 
provide information on magnetic fields integrated along the line of sight. Nevertheless, in conjunction with complementary 
measurements (electron number density), they have be en successfully a pplied to show that the coherence length of magnetic 
fields can be as large as cluster scales, or even bigger I Kim et"al] 1989). S imilarly, magnetic fields have been measured with 
the same amplitudes in high redshift objects at z > 2 ( AthrevZetahf 19981 ) . 

Although magnetic fields are observed on all scales with increasing accuracy, we still do not know where they originated 
from. Primordial magnetic fields are one of the candidates for the origin of magnetic fields in galaxies and galaxy clusters. 
Primordial magnetic fields act as so urces of the primary _B-mode polarisation by generating vorticity in the cosmic plasma 
and additional grav itational waves ( Seshadri fc Subramanian 200 ll : Subramanian et al. 20031 : Mack et al. 2002 ; LewisI 2004 ; 



Tashiro et al. l200d ). However, those waves arise on small scales and therefore do not interfere with the detection of the 



polarisation due to inflation-generated gravitational waves. 

Magnetic fields also create _B-mode polarisation by Faraday rotation, i.e. the rotation of the linear polarisation plane 
through the interaction between CMB photons and magnetised plasma along their path. Faraday rotation distorts the CMB 
polarisation fields like gravitational lensing, so that it induces B-mode polaris ation from the i?-modes. The effects of Fara day 
rotation due to primordial magnetic fields were investigated in some papers i Kosowskv et al. | 2005: Scoccola et al. 2004). In 



addition to primordial magnetic fields, extra-galactic magnetic fields of several /iGauss ( Carilli fc Tavloi l2002h will pr oduce 
Faraday rotation and therefore add to the _B-mode polarisation. This effect was first studied by iTakada et al.l (|200ll ) who 
investigated the _B-mode polarisation due to Faraday rotation in galaxy clusters. They found that homogeneous magnetic 
fields of several /iGauss produce 1 fiK _ B-mode polarisation around I — 1000. However, it has been shown that magnetic fields 
in galaxy clusters are not uniform fsee iMurgia et al.l 12004 and references th erein) . The radial profile of the magnetic fields 



may affect the B-mode polarisation maps on small scales. Ohno et al. ( 20031 ) discuss the possibility of the reconstruction of 



the magnetic fields in a galaxy cluster from the B-mode polarisation map via the Sunyaev-Zel'dovich (S-Z) effect and X-ray 
emission. 

In this paper, we revisit the secondary B-mode polarisation from Faraday rotation in galaxy clusters. We complete the 
study by taking into account galaxies as well. In spiral galaxies, observations show that the magn etic fields spread wid ely not 
only in the galactic plane but also over the halo, with an average amplitude of ab out 10 /xGauss jHumrnel et ah 1991). Even 



in elliptical galaxies, magnetic fields corresponding to many /iGauss are observed ([Greenfield et 



al .1 1 19851 ). 



The paper is organized as follows. In Sec. II, we introduce the formalism of the power spectrum for the Faraday rotation 
angle and we derive the angular power spectrum of B-mode polarisation due to Faraday rotation. In Sec. HI, we discuss the 
two different profiles of gas density and magnetic fields that we use, motivated by observations and numerical simulations. 
In Sec. IV, the resulting B-mode polarisation angular power spectrum is shown, taking into account the effects of both 
galaxies and galaxy clusters. Sec. V is devoted to discussions and summary. Throughout the paper, we use the WMAP values 
of cosmological par ameters, i.e., h = 0.73 (Ho = hx lOOkm/s/Mpc), To = 2.725K, h^Q.^ = 0.0223 and h^Q.^, = 0.128 
i Spergel etahlbood ). 



2 POLARISATION ANGULAR POWER SPECTRUM FROM FARADAY ROTATION 
2.1 Angular power spectrum of the rotation angle 

Our aim is to calculate the angular power spectrum of the CMB polarisation produced by Faraday rotation on the galaxy 
and galaxy cluster scales. To do so, we first need to compute the angular power spectrum of the rotation angle. 
The rotation angle caused by Faraday rotation is given by 

a = „ / ds n^B 7 • b, (1) 
mtuj-^ J 

where u) is the angular frequency, Ue is the free electron number density, B is the magnetic field strength, and 7 and b are 
the directions of the line of sight and the magnetic field, respectively. 

We then u s e the halo-formalism commonly us ed for S-Z power spectrum computations (e.g. Cole fc Kaiser ( 19881 ): 
Makino fc Sutol (|l993h : lKomatsu fc Kitavamal |l999l ll. In this formalism, the angular power spectrum is defined by 



/^ot ^single. 



where Cp^^''^ describes the Poisson term and C/^' 



^single 
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(2) 

gives the contribution from correlated halos. They are expressed as 

(3) 

dn{M, z 



dM 



-b(M,z)ai{M,z) 



(4) 
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where r{z) and V{z) are the comoving distance and the comoving volume, respectively. In Eqs. (|3} and Q, n{M,z) is the 
comoving halo number density of mass M at redshift z and b{M, z) is the linear bias. We adopt the fitting formula given by 
Sheth fc TormenI (|l999l ). 



n{M,z)dM 



M 



1 + 2 



-pr(i/2- 



(1 + 



dv 

V 



(5) 



where v = ^^^|f-jy^ , &c is the critical over density and a is the variance smoothed with a top-hat filter of a scale R = 
{ZM/A-Kp)^/^ and we take p ^ 0.3 a = 0.75 JCoorav fc Shethll2002l ). In Eq. @, we use the linear bias described as b{M, z) — 



l + qv- 1/Sc + 2p/Sc{l + [qvY) jScoccimarro et alJboOll ). 

In Eqs. ([3} and (|4]), ai{M, z) is the projected Fourier transform of the rotation angle obtained in the small angle approx- 
imation, 



ai = 2tv / d69a{e,M,z)Jo{ie), 



(6) 

where oi{9, M, z) is the angular profile of the rotation measurement induced by the magnetic field in a galaxy or a galaxy 
cluster with mass M at redshift z. The rotation angle is obtained from Eq. once a gas distribution and a model for the 
magnetic field in clusters are given. We set 9 = r /Da and 9c = rd Da where Da is the angular diameter distance. To compute 
the power spectrum of the rotation angle, Eq. ([2]), we also need to know the angle between the magnetic field and the line of 
sight. In the following, we assume that the orientation of magnetic fields does not change within a galaxy or a galaxy cluster, 
and we take (|7 • bp) = 1/3. This means that the coherence length of magnetic fields is the virial radius and the orientation 
of the magnetic fields is random from structure to structure. We discuss the case of different coherence lengths in Sec. 14.21 



2.2 B-mode polarisation from Faraday Rotation 



The CMB polarisation can be described using the Stokes parameters (Q, U). If CMB photons travel in the z direction, Q is 
the difference between the intensity in y and x directions, and U is the difference between intensities in directions obtained 
by rotating y and x by 45 degrees. If the CMB radiation passes through a magnetized plasma, it undergoes Faraday rotation 
by an angle a (Eq. H}). The Stokes parameters after passing through the plasma can be written as 



cos 2a sin 2a 
U' ) ^ y -sin 2a cos2q J \ U 

Under the assumption a ^ 1, we get 

Q' + 2aU, 
U' ^U- 2aQ. 

Accordingly, we can write 

Q' ± iU' = (0 -f 2aU) ± i{U - 2qQ) = (1 =F 2ia){Q ± iU). 



(7) 



(8) 



(9) 



The Stokes parameters depend on the choice of the coordinate system. Therefore it is convenient to i ntroduce the rotation 
invariant basis, the -E-modes and the _B-modes (jKamionkowski et al.lll997bl : IZaldarriaga fc Seliakl l 19971 ). These are obtained 
by expanding Q and U in spin-2 spherical harmonics ±2^;"', 



Q±iU = 2^{Ei^ ± iBi^)±2Yr. 

L.m 

The angle of Faraday rotation can also be decomposed as 

l.m 

Applying these decompositions to Eq. ([9]), we obtain dz iB^^ after Faraday rotation as 



(10) 



(11) 



Elm ± iB'im = 



J dfl ±2Yr M =F2i ^ OiimiF™^ j J2 ^Eh_m^±iBl^m^)±2Yi'^^ 
\ l\ ,mi / I2 ,1^2 



(12) 
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From this equation, we can get AS — B' — B, the contribution of the Faraday rotation to the power spectrum of B- 
modes. We give the detailed calculation in the Appendix and we only write the result here. The angular power spectrum of 
the B-mode created by Faraday rotation is written as 

where C^^tfi are the Clebsch-Gordan coefficients, iV, = {2{l - 2)\/{l + 2)!)^/^ and 

K{1, h,l2) = -^ (l' + Li+Ll- 2L1L2 - 2LiL + 2Li - 2L2 - 2L) , (14) 

with L = 1(1 + 1), Li = li{h + 1), and L2 — hih + !)■ Preexisting B-mode polarisation acts like a de-polarisation source 
for the B-modes produced by Faraday rotation (see Appendix, Eq. (j^Zl)- In Eq. (|13|) . we neglect the contributions of such 
preexisting B-modes created, for example, by gravitational waves or gravitational lensing. We checked that these contributions 
are much smaller (two orders magnitude) than primary i5-mode polarisation. 

3 GAS AND MAGNETIC FIELD DISTRIBUTION IN COLLAPSED OBJECTS 

In order to evaluate Faraday rotation on the galaxy or galaxy cluster scales, we need the distributions of the electron number 
density and the magnetic field strength. 

3.1 Gas distribution 

As for the electron number density, we will consider two cases. 



3.1.1 f3 -profile 

We first adopt the /9- profile ( Cavaliere fc Fusco-Femiano 19781 ) for both galaxy clusters and galaxies, 

/ ^2 \ -3/5/2 

ne(r) = ric 1 + — 



(15) 



where r, Vc and Uc are respectively the physical distance from the cluster center, the cluster core radius and the central 
electron density. We explore a different distribution of the electron numb er density in Sec. 13.1.21 Under the assumption of 
self-similarity, Vc oc rvir, ''vir being the virial radius. Cluster observations (|Mohr et al. roughly give rvir ~ lOrc. Using 

the spherical collapse model, the virial radius is related to the mass M of the dark matter halo and the redshift z through 
the following expression, 

1 1/3 



M 



(16) 



where Ac(2) = ISTr^Jlrnz" ''^^ is the spherical over density of the virialized halo ( Nakamura fc Sutol 1997 ). 
In the /3-profile, the number density of baryons can be written as 

/ ^2 \ -3/3/2 

pb{r) = fimpUc ll + ^ j , (17) 

where fi is the mean molecular weight and rUp is the hydrogen mass. We can safely assume that most of the halo mass is 
contained within the virial radius. That is. 



M 



dr Anr ph{T). 



(18) 



From this equation, we can calculate the central electron density 



fib M 



47r 



2^1(3/2, 3/3/2; 5/2; -{r^ir/r.f) 



9.26 X 10"* 



M 



1014 Mq 



^) (^)2Fr^(3/2, 3/3/2; 5/2; -(W'^c)^) 



(19) 



where -Fi(a, /3; 7; z) is the hypergeometric function. Using Vc ~ lOrvir, we find 2^1(3/2, 3/3/2; 5/2; —(rvir/^c)^) is 0.032 for 
P = 0.6 and 0.005 for /3 = 1.0. 
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3.1.2 Navarro-Frenk- White profile 



We also consider a different electron density pro file from the /3-profil e. Namely, we use a density profile motivated by the 
Navarro-Frenk- White (NFW) dark matter profile (|Navarro et al.ll 19971 ') . 
The NFW dark matter density profile is given by: 

Ps 



Pdin(X) = 



(20) 



Here x = r jr^ where is a scale radius, and ps is a scale density at this radius. The scale radius, r-s, is related to the virial 
radius 



c(M,z) 



(21) 



where c is the concentration parameter. In the following, we adopt the concentration parameter of lKomatsu fc SeliakI (|2002l ). 

-0.2 



10 

1 + z 



M 



M,(0) 



(22) 



where M, (0) is a solution to a{M) = Sc at the redshift z = 0. 

T he electron number density profile rie can then be obtained from the NFW dark matter profile following |Komatsu_&_Sel]ak 
( 2002h . For this, three assumptions are made: the gas is in hydrostatic equilibrium in the dark matter potential, the gas den- 
sity follows the dark matter density in the outer parts of the halo, and finally the equation of state of the gas is polytropic, 
Pgas oc p^as (where Pgas, Pgas and 7 are the gas pressure, the gas density and the polytropic index). From these assumptions, 
with Tie oc pgasZ/iTTip, we obtain the electron number density profile 

1/(7-1) 



ne(x) — ric 



1 - F 



1 - 



ln(l + a:) 



where the central electron density ric and the coefficient F are 

-3 



Uc = 3.01 



M 



lO^A/, 



1 Mpc 



\nmJ (1 



(l + c)= 



ln(l + c) 



1+c 



F[c- ln(l + c) 



1/(7-1) 



F = 3r,, 



_i7 - 1 



hi(l + c) 



1 + c 



Komatsu fc SeliakI l|2002h provide the following useful fitting formulae for 7 and 



7 = 1.137 + 8.94 X 10"^ ln(c/5) - 3.68 
7?c = 2.235 + 0.202(c - 5) - 1.16 x 10" 



X 10"'^(c- 5), 



(23) 

(24) 
(25) 

(26) 
(27) 



3.2 Magnetic field distribution 



Magnetic fields in galaxies can reach amplitudes up to 10 pGauss IjHummel et al.lll991f ). Observations bv lMurgia et al.l (120041 ) 
suggest that the distribution of magnetic fields in galaxies and galaxy clusters is such that their strength decreases outward. 
Therefore, in the /3-profile, we adopt the form proposed by those authors: 

-3011/2 



B{r) = Be 1 4- 



(28) 



where Be is the mean magnetic field strength at the center, (3 is the parameter of the /3-profile. From Eq. (|15p . we get B cc n'^^ . 
When /i = 1/2, the energy density of magnetic fields decreases outward like the electron number density, while in the case 
where fi = 2/3, magnetic fields are frozen into matter. 

In the NFW case, the distribution of magnetic fields in galaxies and clusters is given by: 



B{x) ^bAI-F 



ln(l 



m/(7-1) 



(29) 



\i fj, = 2/3, magnetic fields are frozen into the matter whereas, if /i = 1/2, the energy density of magnetic fields decreases 
from the center like the electron number density. 

We further model the growth of magneti c fields in gala xies and galaxy clusters, considering the dynamo process which is 
widely, although not unanimously, accepted (|Widrowll2002l ). In this scenario, magnetic fields are amplified from weak seeds, 
most likely produced by astrophysical processes, by the dynamo process and the time scale of the amplification is of the order 
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Figure 1. The power spectra of Faraday rotation angle caused by 3 /^Gauss magnetic fields in galaxy clusters for the /3-profile case. The 
solid line is for /3 = 0.6 and = 2/3 and the dashed line is for /3 = 0.6 and ii = 1/2. We also plot the case of /3 = 1.0 and = 1/2, and 
f} = 1.0 and /i = 2/3 as a dotted line and a dashed-dotted line, respectively. We use erg = 0.8 and set the CMB frequency to 30 GHz. 



of the dynamical time scale, — r^-^^/GM. Assuming typical magnetic field amplitudes of 3 and 10 /iGauss for clusters 
and galaxies respectively, we model the magnetic field growth by 

^ ^ r 3exp(-(to - t{z))/U) A*Gauss, lO" Mq < M < lO'^'^ M©, 
' \ 10exp(-(to-i(2))Ad) /xGauss, lO" < M < 10^=* Mq, ^ ' 

where to is the present time. 



4 RESULTS 

In order to evaluate the effect of Faraday rotation due to galaxies and galaxy clusters on CMB polarisation, we arbitrarily 
consider that the contribution from galaxies is associated with halos of mass lower than 10^^ Mq. Objects with halo masses 
larger than this limit are considered as galaxy clusters. The particle number density in galaxies is higher than that in galaxy 
clusters. Moreover, in cold dark matter scenario, the number density of galaxies is higher than that of galaxy clusters. Therefore 
the Faraday rotation produced in galaxies presumably cannot be neglected against that generated in galaxy clusters. 



4.1 Rotation angle power spectrum 



In Fig. [TJ we show the angular power spectrum of the Faraday rotation angle due to galaxy clusters only. We assume that 
their mass range is between 10^'^ Mq and 10^^'^ Mq and that all clusters follow the /3-profile with a magnetic field distributed 
following Eq. (I28p and Be = 3 piG. We checked that the Poisson term in Eq. @ dominates in the mass region of galaxy 
clusters. The contribution from correlated halos is thus neglected. We note that increasing (3 increases the amplitude of the 
power spectrum and shifts its peak to smaller angular scales. This is because large l3 value gives steeper profiles and smaller 
core radii. The magnetic field profile also affects the angular power spectrum. However, this dependence is weaker than that 
of the electron density profile. Fig.[5]shows the dependence of the angular power spectrum on t he normalisation of the de nsity 
fiuctuation as- We find l^Ci oc somewhat d i fferen t from the Sunyaev-Zel'dovich (S-Z) case jKomatsu fc Seliakll2002h . We 
compared this result with that of iTakada et al.l (|200ll ) and found it in general agreement. The location of the peaks in our case 
and their case is slightly diffe rent. This can be exp lained by the fact that we take into account distribution of the magnetic 
fields with central peak, while iTakada et al.l l|2001f ) consider homogeneous magnetic fields over the whole cluster scale. 

In Fig. [S] we represent the redshift distribution of Ci, dlnCi/dlnz, for different Is. The power spectrum of the rotation 
angle on large scales (small I) is mostly due to lower redshift galaxy clusters. The rotation angle strongly depends on the 
frequency (see Eq. lU). As a result, low redshift clusters have a larger contribution to C; than in the S-Z effect. For example, 
while most of the contr ibution a,t I — 10000 in the S-Z power spectrum comes from galaxy clusters at redshift z — 2 
i Komatsu fc Seliakll2002l ). ,he power spectrum of the Faraday rotation angle at the same multipole is associated with galaxy 
clusters at redshift z = 1. On scales below / — 100, the contribution to the rotation angle Ci from galaxy clusters at redshift 
lower than z = 0.01 is large. The angular distance of the redshift z — 0.01 is about 40 Mpc while the virial radius of the 
galaxy cluster with mass 10^^ Mq is about 2 Mpc. This means that the small angle approximation used in Eq. © is not 
valid at / < 100. Similarly, we plot the mass distribution of Ci for different Is, dlnCj/dlnM, in Fig. [l] For all Is, the main 
contribution is due to galaxy clusters with masses between 10^^ Mq and lO'^* Mq. The amplitude of the rotation angle power 
spectrum at multipoles larger than I = 5000 is mostly due to clusters with lO'^^ Mq. Comparing the mass contribution of 
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Figure 2. The dependence on the density fluctuation amplitude of the power spectra of Faraday rotation angle. We choose 30 GHz as 
the CMB frequency and adopt the /3-profile with /3 = 0.6, /i = 2/3 and 3 /xGauss magnetic fields. The dashed, the solid, the dotted lines 
are a = 0.9, a = 0.8 and cr = 0.7, respectively. 




Figure 3. The redshift contribution for various I modes. The short-dashed, the solid, the long-dashed and the dashed dotted lines 
represent I = 100, I = 500, I = 1000, I = 5000 and I = 10000, respectively. This figure is calculated for the /3-profile with /3 = 0.6 and 
fi = 2/3 and 3 /iCauss magnetic fields. 



the S-Z effect jKomatsu fc Seliaklliooa '). we find that the Faraday rotation angle is less afTected by the massive halos. This 
is because the high-redshift clusters do not contribute significantly to the amplitude of the Faraday rotation spectrum (see 
discussion on the redshift dependence above) . 



If we take into account the contribution from the whole population of structure including galaxy clusters and galaxies, 
we obtain the results shown in Fig. [S] We calculate both Poisson and halo-halo correlation terms in Eq. @ in the mass range 



1 

13 

a - 

O 0.5 
0.25 



( = 100 
( = 500 - 
( = 1000 - 
I = 5000 

; = 10000 - 



M [Ms 



Figure 4. The mass contribution for various I modes. The short-dashed, the solid, the long-dashed and the dashed dotted lines represent 
I = 100, / = 500, I = 1000, I = 5000 and I = 10000, respectively. This figure is calculated for the /3-profile with /3 = 0.6 and fi = 2/3 and 
3 /iGauss magnetic fields. 
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Figure 5. The angular power spectra of Faraday rotation angle caused by galaxy clusters and galaxies. The solid line is the case of the 
/3-profile with f} = 0.6 and = 2/3 and the dashed-dotted line is the case of the /3-profile with /3 = 0.6 and /i = 1/2. Both are calculated 
using (Tg = 0.8. The dashed and the dotted lines are plotted as the case of erg = 0.9 and erg = 0.7, respectively. In all plots, the CMB 
frequency is 30 GHz. 




Figure 6. The redshift contribution to various I modes. The short dashed the solid, the long dashed and the dashed dotted lines 
represents I = 100, I = 500, I = 1000, I = 5000 and I = 10000, respectively. This figure is calculated for galaxies and galxy clusters with 
the /3-profile for /3 = 0.6 and = 2/3. 



10" Mq < M < 10^^ Mq. In the mass range of clusters lO" Mq < M < lO^*^'^ Mq, we neglect the halo-halo correlation 
term as mentioned previously. The peak of the power spectrum is shifted to smaUer angular scales {I — 20000) because of 
the contribution from galaxies. In Fig. O we use (3 — 0.6 and check that varying fi in Eq. (|28p does not modify significantly 
the power spectrum. We also compare the results with different normalisation parameters. The dependence of the amplitude 
of the rotation measurement is approximately l^Ci oc crf'^^®. The redshift distribution for the power spectrum computed 
using the /3-profile is given in Fig. [S] It also takes into account the evolution of magnetic field Eq. (|3Up which influences the 
redshift distribution. The contribution of high redshift objects (mainly galaxies) is suppressed for z > 1. Moreover the redshift 
distribution at high multipoles is shifted to lower z (see dot-dashed line in Fig. |6}. Most of the Faraday rotation-induced 
signal is thus due to low redshift objects. If we now explore the mass distribution of Ci for different I (Fig. [7|, we notice that 
the contribution from galaxies (left panel) overwhelms that of galaxy clusters (right panel). We also note the difference in 
amplitudes between the two panels is mainly due to the difference in magnetic field strength between galaxy clusters (3 /iG) 
and galaxies (10 jiG) described in Eq. (|30|l . 



Now, following the same approach as for the /3-profile, we compute the angular power spectrum of Faraday rotation 
angle for the NFW profile given in Fig. [8] In this plot, we show the results at 30 GHz, and we assume that the magnetic 
field distribution within clusters or galaxies writes as Eq. (|29|) . For comparison, we plot in the same figure the angular power 
spectrum for the /3-profile and magnetic field given by Eq. (|28|l . The NFW profile is more peaked than the /3-profile. As a 
result, the power spectrum using the NFW profile is shifted to smaller scale (multipoles larger than 30000) as compared to 
those of /3-profile. 
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Figure 7. The mass contribution to various I modes. The left panel shows the contribution from the galaxy mass range while the right 
panel shows that of the galaxy cluster mass range. The short dashed the solid, the long dashed and the dashed dotted lines represent 
I = 100, I = 500, / = 1000, I = 5000 and / = 10000, respectively. This figure is calculated for galaxies and galaxy cluster with the /3-profile 
for p = 0.6 and /i = 2/3. 




Figure 8. The angular power spectra of Faraday rotation angle by galaxy clusters and galaxies at 30 GHz CMB frequency. The solid line 
is the NFW case with magnetic fields frozen into matter, = 2/3, and the dashed line is the NFW case with magnetic fields which energy 
density decrease outward like the electron distribution. We also plot the /3-profile cases with = 2/3 and = 1/2 as the dashed-dotted 
and the dotted lines, respectively. 



4.2 Polarisation power spectrum 

We show, in Fig. (9] the angular power spectrum of the S-mode polarisation caused by Faraday rotation at 30 GHz using 
a /9- profile for the gas distribution. From Eq. (|13|l . the B-mode polarisation by Faraday rotation depends on both the 
primordial i5-mode polarisation Cf and the rotation power spectrum . However, since the rotation power spectrum peaks 
on smaller angular scales than the primordial i5-mode polarisation, the produced _B-mode polarisation reflects the rotation 
angle characteristics. They both peak around I ~ 8000 and for a /3-profile the peak value as a function the main parameters 
of the model is 

We also give the power spectrum of B-mode from Faraday rotation at 100 GHz (thick lines) which is the optimal band for 
CMB observation. 

For comparison, we plot in the same figure the power spectrum of the primary iJ-mode polarisation (thin dotted line), 
the B-mode polarisation due to the gravitational waves (thin dot-dashed line) and the lensing B-mode polarisation (thin 
dot-doted line). The polarisation from gr avitational waves is computed with a tensor to scalar ratio r = 0.6 which is the upper 



limit from the WMAP three years data (jSpergel et al.ll2006i '). This signal dominates on large angular scales. On small scales, 
the power spectrum of the polarisation caused by Faraday rotation dominates the B-mode polarisation due to gravitational 
lensing at frequencies lower than 30 GHz, for reasonable values of the cluster magnetic field. When the frequency increases, 
the gravitational lensing B-mode polarisation dominates. We plot the power spectrum of the Faraday induced polarisation 
for different parameters /3 and ^ describing the electron density and the magnetic field profiles (see line styles in caption). 
Similarly to the rotation angle power spectrum Fig. [1] we note that varying /i has very little effect whereas varying 13 changes 
the amplitude of the power spectrum by two orders of magnitude. 
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Figure 9. The angular power spectra of the B-mode CMB polarisation caused by Faraday rotation for the /3-profile. The solid line is 
for /3 = 0.6 and fi = 2/3, the dashed line is for /3 = 0.6 and fi = 1/2, the dashed-dotted line is /3 = 0.6 and fi = 2/3, and the dotted 
line is for /3 = 1.0 and fj, = 1/2. We use erg = 0.8 and set the CMB frequency to 30 GHz. For comparison, we plot, for 100 GHz, the 
spectrum for /3 = 0.6 and n = 2/3 (thick solid line), and for /3 = 0.6 and fi = 2/3 (thick dashed line). We also give the power spectra of 
the primordial i?-mode polarisation, the B-mode polarisation caused by gravitational waves with a tensor to scalar ratio r = 0.6, and 
the gravitational lensing B-mode polarisation. 




Figure 10. The angular power spectra of the B-mode CMB polarisation by Faraday rotation for the different coherence lengths. The 
solid line corresponds to the model with Ic = rvir- The dashed, the dot-dashed and the dotted lines are for Ic = r'vir/S, Ic = J'vir/4 and 
Ic = Tyir/lO. For comparison, we plot the B-mode polarisation caused by gravitational waves. 



In the above calculation, we assumed that the coherence length of the magnetic field in a galaxy cluster is the virial 
radius. This means that the inner product of the directions of the line of sight and of the magnetic field in Eq. ([ij is 
cons tant over a galaxy cluster. However, m any observations suggest that the coherence length is shorter than the virial radius 
(e.g. [Brandenburg fc SubramanianI l|2005h ). In order to evaluate the eflect of diflerent coherence lengths, we assume that the 



direction of the magnetic field in a galaxy cluster changes smoothly and the inner product of directions in Eq. ([T| has the 
following form, 

7 ■ b = cos (^^^^ I (32) 

where Ic is the coherence length. Taking Ic — r^ir/S, rvir/4 and rvir/lO, we calculate the angular power spectrum of B-mode 
polarisation and plot the results in Fig. 1101 The small coherence length implies many changes of the direction of magnetic 
fields so that depolarisation occurs along the line of sight. When the coherence length is one tenth of the virial radius, the 
amplitude of the angular spectrum is decreased and the peak position shifts to small scales. 

Next, the angular power spectrum of B-polarisation caused by Faraday rotation in galaxy clusters and galaxies is plotted 
in Fig. 111! We plot this figure for a /3-profile with /3 — 0.6 and fi = 2/3, and at a frequency of 30 GHz. The power spectrum 
peaks around I — 10000 with a peak amplitude, which depends on erg, v and Be, 

The _B-mode polarisation induced by galaxies can be the major component of B-mode polarisation on small scales. If the 
average magnetic field of galaxies is 10 /iGauss, the Faraday rotation produces the dominant B-mode polarisation at 30 GHz 
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Figure 11. The angular power spectra of the _B-mode polarisation caused by Faraday rotation in galaxy clusters and galaxies for the 
/3-profile. The solid line is for f) = 0.6 and = 2/3, and the dashed-dotted line is for /3 = 0.6 and ^ = 1/2. We also plot the case of 
(Tg = 0.9 and erg = 0.7 for /3 = 0.6 and /i = 2/3 (dashed and dotted lines respectively). For comparison, we give the i?-mode polarisation 
induced by gravitational waves with r = 0.6. 



for I > 4000. Alternatively, on scales smaller than I = 10000, the _B-mode polarisation generated by galaxies for the same 
magnetic field is expected to dominate at frequencies up to 60 GHz. 

As a galaxy evolves, gas is falling from the halo onto the disk and the halo gas density decreases. The remaining fraction 
of di ffuse baryons in t he galaxy halo at low redshift can be as low as 10 %, whereas in a galaxy cluster it roughly stays equal 



to 1 IjBinnev fc Merrifi eld 19981). The evolution of the gas in a galaxy is expected to affect the power spectrum. We estimate 



the effect of the gas depletion by assuming that only one tenth of the total baryon density is left in the halo and the rest is 
condensated on the galactic disk, so that 



ds fgUeB 7 • b, (34) 



where 



_ / 0.1, M< 1O"M0, 
1, M>1O"M0. 

Ninety percent depletion is the maximum we can expect. Therefore, results in the left panel of Fig. [12] corresponding to 
Eq. (|35|l show the maximum impact of gas depletion on the power spectrum of B-mode polarisation. 

The contribution from galaxies appears only at small scales and is totally overwhelmed by that from galaxy clusters. This 
result tells us that we may neglect Faraday rotation from galaxies, if most of gas is condensed onto the disk. 

In Fig. 111! we showed the B-mode power spectrum taking into account magnetic field evolution accoding to Eq. (|30p . 
However, the magnetic field evolution, especially in galaxy clusters, is not fully understood yet. In some scenarios, the time 



scale of the evolution is shorter than the dynamical time scale (e.g. [Brandenburg fc SubramaniarJ (|2005l )). To highlight the 
importance of the magnetic field evolution, we calculate the angular spectrum in the case of constant mangetic fields and we 
plot the results in the right panel in Fig |12l For higher / modes, the main contribution comes from galaxies at redshifts higher 
than z = 1. Due to this, in the case of constant magnetic fields, the amplitude of the power spectrum is higher and the tail 
at high I is also amplified. 

Finally, we plot the angular power spectrum of the B-mode polarisation induced by Faraday rotation from galaxy clusters 
and galaxies with the NFW profile in Fig. 1131 Here we neglect the gas depletion in galaxies and we consider the magnetic 
field evolution, Eq. (|30p . Using the NFW profile shifts the polarisation spectrum to smaller scales. For /i = 2/3 the spectrum 
peaks at I > 20000, and the value at Z = 10000 is about 



2. ''"-'■'^''^'\w^)\wge-.^ (^6) 



r2 



5 DISCUSSION AND SUMMARY 

In this study, we calculated the secondary B-mode polarisation caused by Faraday rotation. Particularly, we investigated the 
dependence of the B-mode angular power spectrum from Faraday rotation on the electron and the magnetic field profiles in 
galaxies and galaxy clusters. We considered both the /3-profile of electron density as well as an electron density distribution 
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Figure 12. (the left panel) The angular power spectra of the B-mode polarisation caused by Faraday rotation in galaxy clusters and 
galaxies. For the solid line, we assume that 90 % of gas in galaxies are condensed in the disk. The dotted line shows the contribution 
from galaxy clusters only. For comparison, we plot the power spectrum for the case without gas depletion and that due to gravitational 
waves with r = 0.6 as the dashed and the dashed-dotted lines, respectively, (the right panel) The angular power spectra of the B-mode 
polarisation for different evolutions of magnetic fields. The solid line is for the model where magnetic fields are constant. The dotted line 
corresponds to the case where magnetic fields are constant and gas is depleted in galactic halos. We also plot the case where the time 
scale of the magnetic field evolution is the dynamical timescale as the dashed line. The dashed-dotted line is takes into account both 
magnetic field evolution and gas depletion. 
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Figure 13. The angular power spectra of the B-mode polarisation caused by Faraday rotation from galaxy clusters and galaxies at 
30 GHz. The solid line is the case of the NFW profile with /i = 2/3, and the dashed line is for fi = 1/2. For comparison, we plot the 
power spectra for /3-profiles with /3 = 0.6 and fi = 2/3 and /3 = 0.6 and /i = 1/2 (dashed-dotted and dotted lines, respectively). We also 
show the B-mode polarisation induced by gravitational waves with r = 0.6. 

based on the NFW dark matter profile. We modelled the magnetic field structure in galaxies or galaxy clusters motivated by 
observations and further accounted for the redshift evolution of its amplitude. 

We showed that the electron and magnetic field profiles in galaxies and galaxy clusters modify the position of the spectrum 
peak. The amplitude of magnetic fields affects the amplitude of the spectrum, mainly. The B-mode polarisation from Faraday 
rotation at 30 GHz is 0.01 ^K^ at I = 10** for the /3-profiIe, whereas it is 4 x 10~^ /iK^ at the same / for the NFW profile. We 
also investigated the impact of different coherence lengths on the angular power spectrum. Small coherence length induce the 
depolarisation and the peak of the power spectrum shifts to smaller scales and its amplitude is suppressed. 

The Faraday rotation angle power spectrum as well as the B-mode angular spectrum are dominated by the contribution 
from galaxies at redshifts z < 1 and with masses from 10^^ to IO^^Mq. In other words, the detection of a secondary polarisation 
from Faraday rotation would give the average strength of magnetic fields in low redshift galaxies. However, the _B-mode angular 
power spectrum is strongly dependent on the evolution of magnetic fields in galaxies and galaxy clusters, as well as on the 
diffuse gas fraction in galactic halos. If the time scale of the magnetic field evolution is shorter than the dynamical time scale, 
the contribution from galaxies at redshifts higher than z — 1 is dominant. The power spectrum has higher amplitude and the 
peak is on smaller scales. Conversely, if the gas evolution is rapid and most of gas in galaxies condensed onto their disks, the 
contribution from galaxies becomes much smaller and can be neglected. 

The contribution from Faraday rotation from clusters and galaxies appears always small enough, on large scales, to be 
neglected as compared with the primary B-modes from gravitational waves. It is also important to compare the secondary 
B-modes from Faraday rotation with the secondary B-modes caused by gravitational lensing. The B-mode polarisation by 
Faraday rotation dominates that due to gravitational lensing at I > 8000 at 30 GHz. It has been shown that gravitational 
lensing polarisation can be removed from the CMB polarisation map leaving a residual contribution of about only 10% 
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( Hu fc Okamotolbood : iHirata fc Selia^booj ) . Considering this cleaning, the _B-mode polarisation from Faraday rotation will 
exceed the remaining signal from lensing at I = 2000 (30 GHz), or at / = 7000 (100 GHz). These values are obtained in the 
/3-profile case with as ~ 0.8. 

There are other contribution from secondary polarisation sources on small scales that we compare to our results. One of 
them is polarisation from inhomogeneous reionisation. The coupling between the primordial CMB temperature quadrupole 
anisotropy and the density fluctuations of free electrons due to inhomogeneous reio nisation produces a second order polar i- 



sation. Some authors es t imate d this polarisation by using a semi-analytic approach I Liu et al. 2001 : Mortonson fc Hu 20061 ) 
For example. IlIu et al.l l|200ll ) gave l{l + l)Cf /27r ^ 10"^ ^iK^ at / ~ lO". Recentl yjDore et al.l 1I2OO7I) have revisited this 



topic using up-to-date numerical simulation. Their results are consistent with those of Liu et alT i 2001 ). The polarisation from 
Faraday rotation due to the realistic magnetic fields in galaxies and galaxy clusters thus dominates the secondary polarisation 
from inhomogeneous reionisation, at 30 GHz. However, Faraday rotation has a strong dependence on frequency, Ci oc j/^. As 
a result, the polarisation by Faraday rotation becomes subdominant with respect to the polarisation due to inhomogeneous 
reionisation when the frequency exceeds 60 GHz. Next, we considered the polarisation from ionised gas in filaments and large 
scale structure (LSS). T he primary quadr upole produces the polarisation by interaction with the free electron gas captured 
in potential wells of LSS. Liu et aP ( 2005h have evaluated the resulting polarisation and obtained l{l -\- l)Cp /2-k ~ 10~^ /iK^ 
at I ~ 10*. This value is somewhat lower than the polarisation due to Faraday rotation at 30 GHz. While the Faraday rota- 
tion depends on frequency, the polarisation from LSS does not. Therefore, although the polarisation of the Faraday rotation 
dominates the polarisation of LSS at frequencies lower than 30 GHz, the former is subdominant for < 30 GHz. 

Another source of polarisation that needs to be considered is galact ic foregrounds. On small angular scales, the main 
contributions come from synchrotron emission and dust in our Galaxy. iTucci et al.l (|2005j) estimated these contributions, 
focusing on the polarisation at multipoles larger than / = 1000. We extrapolated their results to smaller scales and obtained, 
straightforwardly, /(/ -I- l)Cf /2-k ~ 1 /iK^ at I ~ 10* at 30 GHz. In order to detect the polarisation caused by Faraday 
rotation, we need to remove the galactic foreground polarisation carefully as it will be a major contamination. However, we 
have to bare in mind that such an extrapolation relies on the assumption that there is significant dust contribution on scales, 
/ 10*. On the one hand, present observations indicate that the power spectrum of the dust fluctuations decrease as fast as 

(jMiville-Deschenes et al]|2002h suggesting a fast decrease of the dust contribution on small scales. On the other hand, the 
contribution from synchrotron emission on those small scales depends on the structure and amplitude of the galactic magnetic 
field on the same scales. The distribution of the magnetic fields in the Galaxy has bee n widely studied. Magnetic fields in 
the Galaxy will, therefore, also act as a Faraday rotation source ( Dineen fc Cole j 2005 ). but their contributions need to be 
studied in more detail. 
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APPENDIX A: B-MODE ANGULAR POWER SPECTRUM 

In this appendix, we deal with the derivation of Eq. ([13} from Eq. (|12p . Eq. (|12p involve the integrations of the spin- weighted 
spherical harmonics in the term introduced by Faraday rotation. The integrations of the spin-weighted spherical harmonics 
are calculated through the Clebsch-Gordan coefficients. 



2h + 1 2?2 + 1 ^Im ^1^2 

4tt 21 + 1 



where we use the following property of the spin weighted harmonics. 



l-^il S2-fi2 ^ V 4t V 47!" ^ ^'iml !2m2^il-si !2-S2 Y 2/3 4- 1* '3 ' 



/2/1-fl /2;2 + l ^ ^i,m, w,-., / 47r 

i3"l3S3 

and the orthogonality condition 



,y,7i%Y,7^ = &^^^2Ki2- (A3) 
Substituting Eq. ()AH) to Eq. p2|) . we obtain the representation of B-mode polarisation from Faraday rotation 

-BL = Blm ^ X/ X/ ^ (-l)'l+'2-']y , ^^^^ C'i"^i l2m2<^ifo i22"il™l (-Eiama ~ iBl^^^^) , (A4) 

where we use the following relation to obtain the last equation 

CZ = {-If^'-^CTZ (A5) 

Next, we derive the B-mode polarisation produced by Faraday rotation Ai3;,„. We assume that the rotation fields and 
the primordial polarisation is statistically isotropically and uncorrelated. 
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The angular power spectrum of the B-mode polarisation from Faraday rotation C^^{1) is given by 



(A6) 



(A7) 



where we use the orthogonality of the Clebsch-Gordan coefficients, X/ci,3 '^aZ bfj^aa bp ~ ^cc'S-y-y'- Eq. (|A7|) shows that the 
Faraday rotation produces _B-mode polarisation from iS-mode polarisation and transforms a part of the preexisting B-mode 
polarisation into _E-modes. 

It is hard to calculate Eq. (|A7P directly. To re duce the amount of the c alculations, we rewrite Eq. (|A7P using the following 
equations about the Clebsch-Gordan coefficients (jVarshalovich et al.ll2005l ). 



c: 
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a{a + 1) [c(c -f 1) - a(a + I) + 6(6 + 1)] + 6(6 + I)[c(c + 1) + a{a + 1) - b{b + 1)] o 



2y/a{a + 1)6(6 + l)c(c - l)(c + l)(c + 2) 
After a lengthy computation, we obtain, 

where 

iv, = (2(^-2)!/(^ + 2)!)'/^ 

and 

A-(Z, h,l2) = -^ (i' + i? + Li - 2L1L2 - 2LiL + 2Li - 2L2 - 2L) , 

with L = l(l + 1), Li = li{li + 1), and L2 = ^2(^2 -|- 1). 

To compute the Clebsch-Gordan coefficients in Eq. I|13p . we utilize the approximation of Kosowskv et al. 1 20051 ). 
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This approximation is accurate to only 1% for the worst case a = 6 = c = 2 
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